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The relationship between the applied elastic strain gradient and the induced electric polarization 
in the a-phase polyvinylidene fluoride (PVDF) films under bending conditions has been investigated. 
Our experimental studies have shown that the flexoelectric polarization is linearly proportional to 
the strain gradient and the corresponding direct flexoelectric response is strong. It is reasonable to 
believe that the physical mechanisms behind the flexoelectric effects in polymers and solid dielectrics 
are different. 



The flexoelectric effect, which describes the linear en- 
ergy coupling between the applied strain gradient and the 
induced electric polarization in solid crystalline materi- 
als, was proposed by Kogan in 1964 Kogan's model 
can be written as 



Pi = ^J-ijkl^lSjk 



(1) 



here Pi is the induced polarization; V; is the symbol 
representing the gradient with respect to the direction 
I, fJ-ijki represents the flexoelectric coefficient, and Sjk 
is the applied strain. Initially, Kogan used the term of 
^^piezoelectricity'^ to describe this new phenomenon in his 
paper [ij . Indenbom et al. later suggested using the term 
of ^^flexoelectricity^^ to define the effect since the physical 
mechanisms behind piezoelectricity and flexoelectricity 
are different (H . The term, flexoelectricity^^ , was coined 
by dc Gcnnes for the description of an analogous effect 
in liquid crystals 0. 

The physical picture of the flexoelectricity in solid crys- 
talline dielectrics is clear; it is the inhomogeneous strain 
that breaks the inversion symmetry and induces the elec- 
tric polarization even in centrosymmetric crystalline di- 
electrics 0. The flexoelectric coefficient, fiijki, was pre- 
viously estimated to have an order of magnitude of e/a 
(e and a represent the electron charge and the atomic 
dimension of the unit cell of crystalline dielectrics, re- 
spectively), which is approximately equal to lO^^^C/m 
(see Ref.[l| and references therein). Tagantsev further 
pointed out that this coefficient should be proportional 
to electric susceptibility, which can be written as [H- 



fJ-ijkl = Xijlkli-), 



(2) 



here Xij is the susceptibility of bulk dielectrics under the 
mean field approximation, 7^; is an inhomogeneous sus- 
ceptibility distribution tensor. 

The flexoelectric phenomena in liquid crystals have 
also been extensively studied over the past several 
decades. In 1969, Meyer proposed his phenomenological 
model to describe an analogous flexoelectric effect in liq- 
uid crystals, i.e., the linear energy coupling between the 
splay and bend directing deformations and the induced 
electric polarization, which can be written as Q 



Pd = eii (ndivn) + 633 (n x curln) 



(3) 
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here n is defined as the director that represents the orien- 
tational deformation in liquid crystals, Pd represents the 
induced electric polarization along the direction of n, and 
eii and 633 denote the flexoelectric coefficients related to 
splay and bend distortions in liquid crystals, respectively. 
Like Kogan, Meyer also used the term of "piezoelectric- 
ity^^ in his model when he published his paper. Largely 
based on de Gennes's work 01, we now call such phenom- 
ena in both solid and liquid states as "flexoelectricity^' or 
"flexoelectric effect' . 

The reason that the inhomogeneous deformation and 
the corresponding potential energy could exist in liquid 
crystals and are not destroyed by molecule diffusing in 
the liquid state is that liquid crystals are viscous. It 
is the viscosity that helps liquid crystals preserve the 
aforementioned inhomogeneous deformation, which ren- 
ders the existence of flexoelectricity in liquid crystals. 
Thus the physical mechanism of the flexoelectric phe- 
nomena in liquid crystals is more complicated than that 
in solid crystalline dielectrics. It is well known that 
the flexoelectricity in liquid crystals is closely related to 
the geometrical asymmetry of the mesogen molecules or 
shape polarity , but the relationships between the flexo- 
electric coefficients, en and 633, and molecular structures 
and between the flexoelectricity and the viscosity are far 
from being understood. 
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FIG. 1. (a) a perfect crystal (we use a zig-zag pattern to 
diagrammatically represent the crystalline structure); (b) a 
partially crystallized polymer with both the crystalline struc- 
ture and the amorphous structure. 

The studies of the flexoelectricity in polymers are also 
challenging since the flcxoelectric responses of polymers 
must involve certain time-dependent molecular structure 
changes. The fundamental difference between polymers 
(including solid polymer films) and solid crystalline ma- 
terials resides in the inherent rhcological or viscoelastic 
properties of polymers, i.e., polymers possess both vis- 
cosity and elasticity simultaneously. Therefore, the long- 
term flexoelectric responses of polymers should also pos- 
sess rheological behavior. Besides, the short-term flexo- 
electric responses of polymers are dependent on not only 
their molecule rotation but also the interaction between 
their different phases, which distinguishes the flexoelec- 
tric behavior in polymers from that in both solid crys- 
talline dielectrics and liquid crystals. We use a diagram- 
matic drawing to further demonstrate this viewpoint. In 
Fig. 1, we can see the difference between a perfect crys- 
tal structure (shown in Fig. 1(a)) and a semi-crystalline 
structure (shown in Fig. 1 (b)); a polymer with the lat- 
ter structure can be regarded as a ^^spontaneous compos- 
ite^^ formed with crystalline parts in a glassy matrix i9|, 
which possesses some unique properties. When such a 
semi-crystalline polymer is undergoing deformation, all 
of its parts (both crystalline and amorphous ones) will 
make contribution to its flexoelectric responses. How- 
ever, the net flexoelectricity of the polymer is not the 
sum of that of its different parts or phases; the inter- 
action between those phases might dominate its flexo- 
electric responses under certain circumstances, which has 
been partially confirmed by our recent experimental ob- 
servations [13, Eli- Therefore, the theoretical models and 
experimental techniques developed for the flexoelectricity 
studies of solid crystalline materials and liquid crystals 
might not be appropriate for that of polymers. Unfortu- 
nately, only a few experimental results about the flexo- 
electric properties of polymers have been reported so far; 
our understanding of the flexoelectric behavior in poly- 
mers is still rather limited. 

As shown in Eq. the flcxoelectric coefficient is lin- 
early proportional to the electric susceptibility. Such a 
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FIG. 2. (a) a rectangle shape PVDF film with h — 20mm and 
L = 55mm, the thickness of the film is 13.5/im; (b) the cross- 
section view of an elTective beam system (the rectangle shape 
PVDF film + a cantilever beam); the grounding electrode 
is fabricated on the bottom surface of the PVDF film, and 
then the film is glued by using high strength epoxy to the top 
surface of the cantilever beam; (c) the defiection, w{l), of the 
clamped eff'ective beam system undergoing free vibration, I 
— 11.5mm and L = 55mm; (d) the static defiection, y{l), of 
the clamped effective beam system under an end load F, / — 
11.5mm and L = 55mm. 



relationship has been proposed to exist in polymers too 
Previous experimental studies also showed that 
the order of magnitude of the induced polarization to the 
deformation gradient ratio was approximately lO^^^C/m 
in elastomers [l2j . which is in good agreement with the 
theoretical estimation mentioned, for example, in Ko- 
gan's paper However, Jakli and his co-workers re- 
cently observed a giant flexoelectric phenomenon in cer- 
tain bent-core liquid crystal elastomers (BCLCEs) with 
the nematic phase [lH|, which contradicts the belief that 
the flexoelectric coupling is small in polymers. Prompted 
by their work, we investigated the flexoelectric properties 
of both quasi-amorphous and semi-crystalline polyvinyli- 
denc fluoride (PVDF) films under tensile stretching con- 
ditions; the giant flexoelectric phenomena and other in- 
teresting material behavior have also been observed in 
those films [13, [lli. These experimental observations 
might indicate that the Tagantsev's model given in Eq. 
(12) might not be applied to the flcxoelectric behavior of 
polymers. 

In this Letter, we investigate the flexoelectricity of 
semi-crystalline a-phasc PVDF films under bending con- 
ditions. PVDF is a prototypical semi-crystalline poly- 
mer, which has several forms, such as the apolar a phase, 
the polar /3 phase, etc. The a phase of PVDF has the 
crystallographic point group symmetry 2/m (C2h in the 
Schocnflies notation) , whereas the /3 phase of PVDF has 
the point group symmetry mm2 (C2v hi the Schocnflies 
notation). The method of preparation of the a-phase 
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PVDF film and its XRD cliaracterization have been given 
in Ref. and will not be repeated here. 

The a-phase PVDF films are cut into the pieces with 
the rectangle shape shown in Fig. 2(a). The grounding 
electrode is fabricated on the bottom surface of a rect- 
angle shape PVDF film, and then the film is glued by 
using high strength epoxy to the top surface of a can- 
tilever beam to form an effective beam system shown in 
Fig. 2(b). The direct flexoelectric effect is measured from 



the effective beam system undergoing inhomogeneous de- 
formation generated by its bending movement. Using a 
bending method to measure the flexoelectric behavior of 
solid dielectrics has been extensively studied by Ma and 
Cross (for simplicity, we shall use "Ma-Cross method" 
to represent their bending method from now on) [l6l - [20l | . 
The Ma-Cross method can be summarized as follows: as- 
sume an effective beam system, shown in Fig. 2(c), is 
undergoing a natural or free vibration, its deflection w{l) 
can be written as l2l| 



w{l) = Ai {[sin{l3L) ~ sinh{l3L)] [sin{(il) ~ sinh{pi)] + [cos{pL) + cosh{[iL)] [cos{(il) - cosh{l3l)]) , (4) 



where Ai = C/[sin{l3L) — sinh(l3L)], and we here only 
consider the fundamental mode of the vibration. There- 
fore, we have /3L = 1.875 [2l[. C can be determined from 
the boundary condition, i.e., the measured displacement 
in the Z axis direction of the effective beam system. Both 
the direction of w{l) and the Z axis are shown in Fig. 
2(c). 
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FIG. 3. (a) a lock-in detection setup; C: current measurement, 
FG: function generator, MS: MicroStrain DVRT probing sen- 
sor, DLT: dynamic loading tester; (b) the top electrode is 
fabricated on the top smface of the PVDF film; I — 11.5mm 
and L — 55mm. 

If we further assume that the deflection of the effective 
beam system under bending conditions is infinitesimal, 
then the relationship between its strain in the X axis 
direction and its defiection 'w{l) in the Z axis direction 
can be written as 1221 



^11 



(5) 




34.5 ± 2.0 nC/m 



Strain Gradient (m ') 

FIG. 4. The effective flexoelectric coefficients, Hcb, of a-phase 
PVDF films measured at room temperature; the line with the 
"circle maker", whose slope is ficH ~ 34.5 ± 2.0/iC/m, repre- 
sents the measured results by using our method; the line with 
the "square maker", whose slope is ficit = 37.3 ± 2.0/iC/m, 
represents the measured results by using the Ma-Cross's 
method. 



where Sn is the applied strain; "1" represents the X axis, 
I is the variable on the X axis. Similarly, we also use "2" 
and "3" to represent the Y axis and the Z axis, respec- 
tively. The minus sign given in Eq. ([5]) only indicates 
whether the top surface of the beam system is in the ten- 
sile state (^ii > 0) or the compressive state {Sn < 0). 
The absolute value of the gradient of the strain along the 
Z axis direction is then given by 



9^11 



dz 



(6) 



According to the Kogan's model defined in Eq. ([T]), 
the linear relationship between the elastic strain gradient 
and the electric polarization of the PVDF film under the 
bending condition, shown in Fig. 2(c), can be written as 



Ps ^ f^eff 



dSn 
dz 



(7) 



where P3 represents the induced electric polarization 
along the Z axis direction; fJ-eff is defined as the effective 
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flexoelectric coefficient . P3 can be further calculated by 
using the following equation: 



27rfA' 



(8) 



where / is the induced alternating current; / is the driv- 
ing frequency of a dynamic loading tester (DLT) dia- 
grammatically shown in Fig. 3(a), and A represents the 
area of the top electrode, which is demonstrated in Fig. 
3(b). In our experiment, the film samples we studied were 
neither stretched nor poled so that we could exclude pos- 
sible extrinsic factors from our flexoelectric measurement; 
besides, such samples can be approximately assumed to 
possess a point group symmetry cxDoom since the a-phase 
domains in those samples are randomly oriented. Thus 
we can rc-definc fJ-eff as: /ie// = J^Mii + (1 + 1^)^12, here 
v represents the Poisson's ratio, and both /xn and mi 2 
are the flexoelectric coefficients in matrix notation [23|. 
Since determining /in and /ii2 needs tedious measure- 
ment of the applied strain gradient and the induced elec- 
tric polarization and their explanations also need more 
space, which is not suitable for a brief research report, 
we will not discuss the details of the measurement here. 
In this study, we mainly investigate the effective flexo- 
electric responses of a-phase PVDF films under bending 
conditions. 

There is still a concern about the Ma-Cross method, 
i.e., whether the assumption that the effective beam sys- 
tem is undergoing a natural vibration is appropriate or 
not. This is because Ma and Cross actually used the 
forced vibration not the natural vibration to generate 
the bending movement of the cantilever beam with one 
end clamped in their studies [l6l-|20j. Therefore, the de- 
flection w{l) defined by Eq. ^ should be re-examined, 
which is another research objective of this Letter. We 
consider a clamped effective beam system under an end 
load shown in Fig. 2(d); its static deflection, y{l), is given 
by HI 



y(i) 



(3L 







(9) 



where A2 = —F/6EI; F represents the end load, E and / 
are the modulus of elasticity and the moment of inertia of 
the clamped effective beam system, respectively; I and L 
are defined in Fig. 2(d). If we let F vary slowly, y{l) will 
then become quasi-static, which means that the bend- 
ing movement of the beam system shown in Fig. 2(d) 
could be generated by a forced low-frequency vibration. 
Similarly, A2 can also be determined from the boundary 
condition, i.e., the measured displacement in the Z axis 
direction of the effective beam system. Both the direc- 
tion of y{l) and the Z axis are shown in Fig. 2(d). When 
y{l) is very small, the relationship between the gradient 
of the generated strain along the Z axis direction and the 
deflection y{l) can be written as 



dz 



(10) 



Our flexoelectric measurement is summarized as fol- 
lows: the aforementioned forced low-frequency vibration 
of the beam system is driven by the DLT shown in Fig. 
3(a), which is modified from a commercial piezo actua- 
tor; a function generator (FG) is used to provide an AC 
sinusoidal signal with a frequency of / = O.SHz to an am- 
plifier to control the periodic movement of the DLT and 
to a lock-in amplifier (LIA) as the reference signal. The 
dynamic bending of the beam system is controlled via 
the DLT, and its displacement or deflection in the Z axis 
direction is detected by using a MicroStrain DVRT prob- 
ing sensor; the measured displacement values are used to 
calculate the strain gradient along the Z axis direction 
via Eqs. (jU and ([6|) (by using the Ma-Cross method) or 
via Eqs. (jH]) and ([TU]) (by using our method). Finally, 
using Eq. ([5]), we can obtain the induced polarization, 
P3; here / is collected by using the LIA; I, L, and A are 
shown in Fig. 3(b). Under the same physical conditions 
and at room temperature, we measured the flexoelectric 
responses of a-phase PVDF films under bending condi- 
tions. We obtained approximate values of the effective 
fiexoelectric coefficients, fieff = 34.5 ± 2.0/xC/m (via our 
method) and /ios = 37.3 ± 2.0/xC/m (via the Ma-Cross's 
method). The linear relationships between the induced 
polarization values and the applied strain gradient values 
are shown in Fig. 4. 

We have to emphasize that, in both the Ma-Cross's 
method and our method, it is critical to keep the gen- 
erated deflection as small as possible since Eq. ([5]) and 
Eq. (fTO)) can only be used to calculate the strain gradi- 
ent when the deflection is infinitesimal. So, in theory, 
the testing point, the spot covered by the top electrode 
shown in Fig. 3(b), on the beam system should be chosen 
near its clamped end. However, a perfectly rigid clamped 
condition for a cantilever beam can never be achieved 
in practice; an imperfect clamped condition often makes 
the lock-in detection very difficult. Therefore, the test- 
ing point for such flexoelectricity measurement should be 
carefully chosen to reduce the inevitable errors caused 
by the imperfect clamped condition. In our experiment, 
we chose the spot on the beam system, where the ratio 
l/L w 1/5 shown in Fig. 3(b), as our testing point. 

In conclusion, a giant direct flexoelectric response 
in a-phase PVDF films under bending conditions has 
been observed; we also found that the induced fiexo- 
electric polarization is indeed linearly proportional to 
the applied strain gradient in our samples. It is per- 
haps interesting to compare the giant flexoelectric phe- 
nomena observed in polymers with that in solid di- 
electrics. Ma and Cross observed a giant flexoelectric 
phenomenon, fii2 ~ 100/xC/m, in a barium strontium ti- 
tanate (Bao.eySro.asTiOa or EST) composition [l^; in 
their studies, the EST composition is in the paraelectric 
state just above its Curie temperature, and, therefore, 
its relative permittivity value is very large, £r ~ 16000 
p^ . According to the Tagantsev's model, it is not sur- 
prised that giant flexoelectricity can be observed in such 
a solid material. However, the average relative permittiv- 
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ity value of our a-phasc PVDF samples is only Sr ~ 11, 
which is three orders of magnitude smaller than that of 
the BST. This indicates that the Tagantsevs model might 
not be applied to the flexoelectric behavior of polymers. 
Finally, we also verified that the measured fiexoelectricity 
results by using the Ma-Cross method and our method 
are in good agreement when the frequency of the periodic 



external driving force is low. 
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